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Abstract 



We study in this paper chameleon cosmology applied to Friedmann-Robertson-walker space, 
which gives rise to the equation of state (EoS) parameter larger than —1 in the past and less than 
— 1 today, satisfying current observations. We also study cosmological constraints on the model 
using the time evolution of the cosmological redshift of distant sources which directly probes the 
expansion history of the universe. Due to the evolution of the universes expansion rate, the model 
independent Cosmological Redshift Drift (CRD)test is expected to experience a small, systematic 
drift as a function of time. The model is supported by the observational data obtained from the 
test. 

PACS numbers: 04.50.h; 04.50.Kd 

Keywords: Chameleon cosmology; Phantom crossing; Bouncing universe; Cosmological Redshift Drift. 



* Electronic address: |hosseinf@guilan.ac.ir| 
^Electronic address: a. salehi@guilan. ac . ir 



1 



1. INTRODUCTION 



The recent astrophysical data indicate that there is a dark energy (DE) providing approx- 
imately two third of the current universe energy density which explain the current cosmic 
acceleration l|j2|- The most obvious candidate to explain DE is the cosmological constant 
which can fit observations well. However it is so small ( of order 10~^^eV), in comparison 
with the Planck scale {10^^GeV7) that suffer from fine-tuning and the coincidence problems 
{sl- Numerous other DE models are produced by some exotic matter like phantom (field 
with negative energy) or some other (usually scalar ) matter ^. Unfortunately, such scalar 
field is usually very light and its coupling to matter should be tuned to extremely small val- 
ues in order to be compatible with the Equivalence Principle. In a sense, the cosmological 
evolution of the scalar field contradicts with the solar system tests 5|. 

From particle physics point of view, there is also wide-spread interest in the possibility 
that, in addition to the matter described by the standard model, our universe may be 
populated by one or more scalar fields. This might be a feature in high energy physics 
beyond the standard model and are often related to the presence of extra-dimensions j^. 
The existence of scalar field may also explain the early and late time acceleration of the 
universe It is most often the case that such fields interact with matter; directly due 

to a Lagrangian coupling , indirectly through a coupling to the Ricci scalar or as the result 
of ,ua„t.n .OOP co^ect,™. Q-Q. If tKe scala. field .e.f-.nt.action. are negligible, then 
the experimental bounds on such a field are very strong; requi ring it to either couple to 
matter much more weakly than gravity does, or to be very heavy 17|- 20|. Recently, a novel 
scenario was presented that employed self-interactions of the scalar-field to avoid the most 



restrictive of the current bounds 



211 ]. In the proposed model, a scalar field couples to matter 



with gravitational strength, in harmony with general expectations from string theory whilst 



at the same time remaining very light on cosmological 



scales. In this paper we will go further 



and, contrary to most expectations presented in |2l(], allow scalar field which is very light 
on cosmological scales, to couple to matter much more strongly than gravity does, and yet 
still satisfies the current experimental and observational constraints. The cosmological value 
of such a field evolves over Hubble time-scales and may cause the late-time acceleration of 



our universe 



22I ]. The crucial feature these models possess is the mass dependency of their 



scalar field on the local background matter density. On the earth where the density is about 



2 



10 times higher than the cosmological background, the Compton wavelength of the field 



is small enough to satisfy all existing tests of gravity 



On the other hand, in the solar 



system, where the density is several orders of magnitude smaller, the Compton wavelength 



of the field can be much larger 



24|. This means that, in these models, the scalar field may 



have a mass in the solar system much smaller than was allowed. In the cosmological scale, 
the field is lighter with its energy density evolves slowly over cosmological time-scales and 
it may be considered as an effective cosmological constant. Although the idea of a density- 



dependent mass term is not new 25|]-30|, the work presented in 21 1 



22l | is outstanding in 



that the scalar field can couple directly to matter with gravitational strength. 

There are cosmological models 3l| were initially proposed to obtain a model of dark 
energy with the EoS parameter u > —1 in the past and < — 1 at present. These models 



can be viewed as dynamical model for dark energy with the feature that t 



heir 



can smoothly cross over the cosmological constant barrier u = —1 32j 33|. To construct 



]oS parameters 



the model it is necessary to add extra degrees of freedom with un-conventional features to 
the conventional single field theory if we expect to realize viable models in the framework of 
gravity theory 3^. We believe that the chameleon model investigated in this work is capable 
to do the job. It predicts u crossing scenario in the history of the universe in addition to 
the bouncing behaviour of the universe. 

Since the cosmic acceleration affects the expansion history of the universe, to understand 
the true nature of the driving force, mapping of the cosmic expansion of the universe is 
very crucial 35|. Hence, we require various observational probes in different redshift ranges 
to understand the expansion history of the universe. The observational tools for probing 
the cosmic acceleration broadly fall into two categories. Geometrical probe that deals with 
large scale distances and volume which include luminosity distance measurements of SNe la, 
angular diameter distance from first CMB acoustic peak, baryon acoustic oscillations and 
so on. Dynamical probe that investigates the growth of matter density perturbations to 
give rise to the large scale structure such as galaxies and clusters of galaxies in the universe. 
Using supernova as standard candles is a popular method of constraining the properties 
of dark energy. Although these methods are very simple and useful in constraining the 
various dark energy models, at present, the luminosity distance measurements suffer from 
many systematical uncertainties such as extinction by dust and gravitational lensing 



Measuring the expansion history from growth of matter perturbations also has its limitations. 
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It requires prior information of exact value of matter density, initial conditions, cosmological 
model and so on [s^js^]. So the question arises, in any probe which is simple, depends on 
fewer priors and assumptions. 

One possible probe is Cosmological Redshift Drift (CRD) test which maps the expansion 
of the universe directly. It assumes that the universe is homogeneous and isotropic at the 



cosmological scales 38] . The CRD test while is based on very simple and straightforward 



physics, observationally it is a very challenging task and requires technological breakthroughs 



39|. This probe measures the dynamics of the universe directly via the Hubble expansion 



factor. The time drift of the cosmological redshift probes the universe in the redshift up 
to 2; = 5, whereas the other cosmological tests based on SNela, BAO, weak lensing or 
number counts of clusters have not advanced beyond z = 2. Its other advantage is that it 
has controlled systematical uncertainties and evolutionary effects of the sources 40|. The 



test was first proposed by the author in 4l|| in terms of the scale factor. So, the universe 
expansion rate is directly measured by the time evolution of the scale factor, or change 
in redshift, z. The measurement of the redshift of an object today, is different from its 
measured value after a time interval of several years. As explained in 42[ the redshift drift 



signal, i;, is very small. In 



43| . author was the first to suggest the possibility of measuring 



the redshift drift by observing Lyman a {LYa) absorption lines in the spectra of quasars 
(QSOs), which reinforcs the importance of this probe. 

This work is arranged as follows. In the next section we present the chameleon model. 
We also investigate the conditions for the EoS parameter in the model to across —1 and 
moreover examine the possible bouncing of the universe. In section three, by using the field 
equations we introduce the observational data and describe their inclusion in our analysis 
through driving CRD test for the model. We also revisit the Chevallier-Polarski-Linder 
(CPL) model in a comparison with our model. In Section four, we present the summary 
and remarks. 
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2. THE MODEL, PHANTOM CROSSING AND BOUNCING UNIVERSE 



We consider the chameleon gravity in the presence of cold dark matter with the action 
given by, 

S = J - ^0,^0'^ + V{<P) + /(0)£,„]v^dx^ (1) 

where R is Ricci scalar, G is the newtonian constant gravity and is the chameleon scalar 
field with the potential V{(f)). Unlike the usual Einstein-Hilbert action, the matter La- 
grangian Cm is modified as f{(f))Cm, where /(0) is an analytic function of 0. The last term 
in the Lagrangian brings about the nonminimal interaction between the cold dark matter and 
chameleon field. The variation of action ([1]) with respect to the metric tensor components 
in a spatially fiat FRW cosmology yields the field equations, 

3H' = Pmf + ^^' + V{^), (2) 

2H + 3H^ = -jp„J - ^0' + ^(0), (3) 

where we put SttG = c = h = 1 and assume a perfect fiuid with pm = 'jpm ■ The energy 
density pm stands for the contribution from the cold dark matter to the energy density. Also 
variation of the action ([T]) with respect to the scalar field provides the wave equation for 
chameleon field as 

4> + 3H<p = -V' -jpmf', (4) 

where prime indicated differentiation with respect to 0. From equations (E]) and @, 
one can easily arrive at the relation (extended conservation equation) 

[pmf) + 'iHpm{l + l)f = lpm(t)f\ (5) 

which readily integrates to yield 

^"^ /(l-7)a3{l+7)' 

where M is a constant of integration. From equations (j2]) and ([3]) and in comparison with 
the standard friedmann equations we identify pe// and pe// as 

Pe//=P,n/+^0' + V^(0), (7) 



1 



Peff =7Pmf + 7^(1) -V{(f)), 



(8) 



with an effective equation of state, Pe// = ^ejfPefj- From equations ([7]) and®, we yield 

2 



eff 



-^[Pmfil - ^eff) - V{uJeff + 1)]. 



(9) 



The possibility of phantom crossing can be discussed in here. If crossing happens at time 
t = tcrossi from equation (9) we have 



-(7 + l)p^/. 



(10) 



Equation (11 01) shows that, for phantom crossing, / must be negative. Also equations ([7]) , 
^ and © show that, if (fP' > —(1 + 7)pm/ then u^jj < —1 and if 0^ < —(1 + 7)pm/ then 

Ueff > -1. 

Now in order to close the system of equations we make the following ansatz: We take 
the potential to be power law in as V{(j)) = Vocf)^ and assume / behaves exponentially as 
/(0) = /qc'"^ where n, 6, Vq and /o are arbitrary constants. Figure (la) shows that the possi- 
ble crossing happens only when / is always negative. In order to keep / negative a runaway 
behavior is considered for it which also satisfies equation f llOp . In figure (lb), 0^ never 
intersect — (l+7)pm/ as / selected to be always positive so phantom crossing does not occur. 



-X a) 




\ \ 




-(1+Y)p > ^^^^ 
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eff 








Fig.l: Plots of We//, 0^ and -(1 + 7)p^/ with f = fo e^-^and V = V"o0",(& = -1 ,n 
for a) / < with {fo = -7) and b) / > with (/o = 2). 
Initial values are 0(0) = 1, 0(0) = -0.8, a(0) = 0.8 and d(0) = 0.1. 
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Alternatively, we may take / to be power law in 0, = fo4>"', and potential behaves 
exponentially, V{(j)) = Vqc*'^. There are no priori physical motivation for these choices, so it 
is only purely phenomenological which leads to the desired behavior of phantom crossing. 

In figure (2a), we examine phantom crossing for exponential / and power law potential. 
We see that, as expected for / < we obtain crossing and for / > do not. At / = that 
0^ vanishes, from equation ( ITOj) and figure (2), Ueff just becomes tangent to the hne — 1. 
The same argument can be made in figure (2b) for power law / and exponential potential. 



Fig.2: Plot of Ueff with a) f = fo e^'^and V = VqcP" or b) ^ = Vg e^<^and / = /o0" 
for fixed Vq and three different value of /o = —7 ,/o = and /q = 2,(6 = — 1 ,n = —1). 
Initial values are 0(0) = 1, 0(0) = -0.8, a(0) = 0.8 and a(0) = 0.1. 

By choosing t ^ to be the bouncing point, the solution for a{t) and H{t), shown in figure 
(3), provides a dynamical universe with contraction for t < —0.02, bouncing at t = —0.02 
and then expansion for t > —0.02. 



a) 
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Fig. 3: Plots of scalar factor a{t) and Hubble parameter H, at the bounce for 
fg Q^'t'^ /o = -7 and h = -1. Initial values are 0(0) = 1, 0(0) = -0.8, a(0) = 0.8 and a(0) = 0.1 

A detailed examination on the necessary conditions requires for a successful bounce shows 
that during the contracting phase, the scale factor a{t) is decreasing, i.e., d < 0, and in the 
expanding phase we have d > 0. At the bouncing point, d = 0, and so around this point 
d > for a period of time. Equivalently in the bouncing cosmology the Hubble parameter 
H runs across zero from i7 < to iJ > and if = at the bouncing point. A successful 
bounce requires that the following condition should be satisfied around bouncing point. 

From figure (2a) and (2b) for / < 0, we see that at t — )■ 0, Wg// becomes less than —1 and 
H is positive which satisfies the above condition. Also we see that at the bouncing point 
where the scale factor a{t) is not zero we avoid singularity faced in the standard cosmology. 

3. COSMOLOGICAL REDSHIFT-DRIFT TEST(CRD) 

The observed redshift of a distant source is given by 

where tg is the time at which the source emitted the radiation and tobs is the time of ob- 
servation. In expression f|T2|) . any peculiar motion of the object is ignored. After the time 
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interval Atobs the source redshift becomes 

.(u,+A...)^ °'^'-;fy -i. (13) 

where Atg is time interval of the source emitted. One can write the first order approximation 
of the above equation as 

Az (aitobs) - a(ts)) 



or in term of the Hubble parameter: 



(14) 



Ho[l + z-^]. (15) 



Equation ( !T5|) is also known as McVittie equation 4J]. It clearly shows that z traces H(z), 
the Hubble parameter at redshift z. As stated in the introduction, z measures the rate of 
expansion of the universe: i > (< 0) indicates the accelerated (decelerated) expansion of 
the universe, respectively. For a coasting universe z = 0. The redshift variation is related 
to the apparent velocity shift of the source: 

Az 

Av = c^=—. (16) 

l + z ^ ^ 

Thus, using equation (fT5|) and (fT6|) one can obtain, 

,17) 

where i) = and Hq = lOOhKm/ sec/Mpc. In the standard cosmological model 

[AC DM), the change in redshift for a time interval of Atobs = lOyr is Az ~ 10~^. For 
a source at redshift z = 3, the correspondin g ve locity shift is of the order of Av ~ 7.5cm/ s. 



To measure this weak signal, the author in 



43| pointed out the detection of signal of such 



a tiny magnitude might be possible by observation of the LYa forest in the QSO spectrum 
for a decade j4o| . 

To observe such a tiny signals, a new generation of Extremely Large Telescope (ELT), 
equipped with a high resolution, extremely stable and ultra high precision spectrograph 
is needed. Using the Cosmic Dynamics Experiment (CODEX) operation and performing 



Monte Carlo simulations of quasars absorption spectra 45 
ments. 



one obtains the z measure- 
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In the following, we use three sets of data (8 points) for redshift drift generated by 
performed Monte Carlo. We investigate various models of the universe which can explain 
late time acceleration, and compare the chameleon scalar field model with these models and 



observational data 



45|-|47| 



3.1. ACDM model 

One of the widely studied model of dark energy is XCDM parametrization in which the 
dark energy is characterized by time independent EoS parameter, with p = u^p. Here 
ACDM model means standard cosmological model with u = —1, ^2^0 = 0.3 and Qao = 0.7. 
For acceleration in this model we need a;^,. < The Friedmann equation in this model for 
a fiat universe is given by 



Ha 



(18) 



where and are the fractional matter and dark energy densities at the present epoch 
respectively with + fix = ^■ 

3.2. CPL model 



Another popular parametrization which explains evolution of dark energy is the CPL 
model 49[ in which in a fiat universe the time varying EoS parameter is parametrized by. 



uj{z) = Uo + UJli— — ). 

1 + z 

The Hubble parameter in the model is given by, 

"^^"^^^ ' = + zf + (1 - nj{l + z)3(l+-o+-i) X 



(19) 



exp 



1 + z 



(20) 



The parametrization is fitted for different values of Uq and Ui. The velocity drift with 
respect to the source redshift is shown in figure (4) [48|. As can be seen, the best fit values 
are for uq = —2.4 and Ui = 3.4. 



10 




3.3. Chameleon model 

By using a = and taking derivative we get 

^ = + , )^,(,)^. (21) 
From equations and (^T^ we obtain the following expression for EoS parameter, 

^e// = -l+ ^ 3/ , (22) 

where r = ^ and r^*^^ = In figure (5), it can be seen that for different values of Vq, the 
EoS parameter crosses —1 at different z. It shows that for larger values of Vq, the phantom 



crossing is more appropriate with the observational data 50 1. 
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Fig.5: Plot of uJeff{z) with f = fo e^-^^^^and V = for fixed value /o 

{b = — 1 ,n = —1). Initial values are 0(0) = 1, 0(0) = —0.8. 



-7. 



It also shows that Wg// will decrease before crossing along with the increase of redshift 
z. This result implies that the behavior of dark energy is different in different slices of 
redshifts, and inspires us to separate redshifts into several pieces and to investigate each 
piecewise separately which can be investigated in another work. 

From equation (122|) we also obtain, 

■H{z)^^ r„ r 1 



exp 



eff\ 



-dz 



(23) 



'0 1 + 5 

From numerical computation of uJeff{z), one can obtain H{z). Then, using equation ([1] 
one leads to the velocity drift. Figure (6) shows the velocity drift against z for various 
potential functions. In comparison with the CPL model our model is in a better agreement 
with the experimental data for z > 3 whereas the CPL model is more appropriate for z < 3. 
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4. CONCLUSION 

In this paper, we study the evolution of the gravitational and scalar fields in chameleon 
cosmological model in which a light scalar field (chameleon field) nonminimally coupled to 
the matter Lagrangian. We find that the evolution of the scale factor of the universe is 
non-singular in a bouncing cosmology, with an initial contracting phase which lasts until 
to a non-vanishing minimal radius is reached and then smoothly transits into an expanding 
phase which provides a possible solution to the singularity problem of standard Big Bang 
cosmology. In addition. The evolution of the cosmological EoS parameter, with a transition 
from a; > — 1 in the past to the uj < —1 in the recent past is favored for negative /(</>), in 
agreement with the current observational data. The phantom crossing occurs for different 
values of Vq within the range of observationally accepted redshift z (0.1 < z < 0.45). 

We then analyze the chameleon model with the CRD test. The variation of velocity drift 
for different values of Vq with redshift z is shown in figure (6). A comparison between this 
model with CPL model shows that it is in better agreement with the experimental data for 
values of redshift z greater than 3. We also obtain that the Hubble parameter, H(z), in the 
model is potential dependent across the redshift. In order to better constrain the parameters 
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of the model we need more redshift drift data in the redshift range z < 1.9. 
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